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ABSTRACT 



We present simultaneous high-temporal and high-spectral resolution observations at optical and soft X-ray wavelengths of the nearby 
flare star CN Leo. During our observing campaign a major flare occurred, raising the star's instantaneous energy output by almost 
three orders of magnitude. The flare shows the often observed impulsive behavior, with a rapid rise and slow decay in the optical and 
a broad soft X-ray maximum about 200 seconds after the optical flare peak. However, in addition to this usually encountered flare 
phenomenology we find an extremely short (j iec as 2 sec) soft X-ray peak, which is very likely of thermal, rather than non-thermal 
nature and temporally coincides with the optical flare peak. While at hard X-ray energies non-thermal bursts are routinely observed 
on the Sun at flare onset, thermal soft X-ray bursts on time scales of seconds have never been observed in a solar nor stellar context. 
Time-dependent, one-dimensional hydrodynamic modeling of this event requires an extremely short energy deposition time scale Td ep 
of a few seconds to reconcile theory with observations, thus suggesting that we are witnessing the results of a coronal explosion on 
CN Leo. Thus the flare on CN Leo provides the opportunity to observationally study the physics of the long-sought "micro-flares" 
thought to be responsible for coronal heating. 
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1. Introduction 

The basic energy release in solar and stellar flares is thought 
to occur in the coronal regions of the underlying star, by con- 
version of non-potential magnetic e nergy through magnetic re- 
connection dPriest and F orbes, 2002). Particle acceleration takes 
place and the accelerated particles and/or plasma waves move 
along the mag netic field lines and reach and heat coole r atmo- 
spheric layers (ISvratovskii and Shmeleval Il972t iBrownl Il973l) . 
Because of their small cooling times the heated photospheric 
layers start radiating immediately, serving as a proxy indicator 
for non-thermal particles, while the heated dense chromospheric 
layers "evaporate" leading to a soft X-ray flare dBrownl Il973t 
lNeuperdfl97iaiPerei l2~000). The heating and cooling time scales 
and the dynamic response of the various atmospheric layers are 
important for our u nderstanding of solar a nd ste l lar flares and 
of coronal heating dKostiuk and Pikelneti Il974t ISomov et all 
I198U iTamres et all 119861: iFisheti 11981) : a popular hypothesis 
dKlimchukL 120061 iParkeil Q988) even attributes all of the re- 
quired energy input to the quiescent solar corona to (nano)-flare 
input. 

Solar and stellar flares are observed to occur over 
vastly different time scales ranging f rom a few seconds 
dVilmer et all 1 1994 [Schmitt et all ll 993b to more than a week 
dKiirster and Schrnitn, Il996b . and the response of the heated 
plasma sensitively depends on the mode of energy input. Since 
the flare process involves plasma differing in temperature and 
density by more than four orders of magnitude, varying on short 
spatial and temporal scales, simultaneous multiwavelength data 
with sufficient spectral coverage and sufficient spectral and tem- 
poral resolution are required for any sensible observational di- 
agnostics. Since the primary energy release and energy loss 



Correspondence to: J.H.M.M. Schmitt, jschmitt@hs.uni-hamburg.de 



process of a flare is coronal, space-based X-ray observations 
are required. Using ESA's XMM-Newton X-ray observatory 
and ESO's Ultraviolet- Visual Echelle Spectrograph (UVES) on 
Kueyen on May 19'' 1 2006 we observed a spectacular flare on 
the nearby (d = 2.39 pc) flare star CN Leo (spectral type: M5.5- 
M6.0, T ef f = 2800-K 2900 K). Desp ite its low apparent rota - 
tional velocity of v sin i < 3.0 km s -1 (Fuhrmeis ter et all l2004). 
CN Leo shows all the attributes of a magnetically active star, in- 
cluding chromospheric and coronal emission as well as flarin g 
in the optical and in the X-ray bands (Fuhrmeist er et alll2007b . 

2. Observations and data reduction 

The XMM-Newton observatory carries three co-aligned X-ray 
telescopes and a co-aligned smaller optical telescope equipped 
with the Optical Monitor (OM); this satellite and its instru- 
ments are described in detail in a Special Issue of Astronomy 
& Astrophysics Vol. 365 (Jan. 2001). For the CN Leo observa- 
tions we used the medium filter in Full Frame and Large Window 
mode for all EPIC observations. The OM was operated in fast 
window mode, allowing U band flux measurements with a time 
resolution of 1 sec, limited by the chosen read-out sequence. 
The X-ray data were reduced with the XMM-Newton Science 
Analysis System (SAS) software, version 7.0. The basic XMM- 
Newton data consist of individual photons with known arrival 
time, position and energy; using standard filtering criteria these 
data are further processed into EPIC lightcurves and spectra. 

UVES is an cross-dispersed echelle spectrograph mounted 
on the Nasmyth B focus of the Kueyen telescope (UT2) cover- 
ing the wavelength range from about 3000 A to 1 1 000 A with a 
typical spectral resolution of ~ 40000. The temporal resolution 
of the high-resolution spectra is determined by the sum of ex- 
posure and read-out times (a few minutes in our case), a small 



2 J.H.M.M. Schmitt, F. Reale, C. Liefke, U. Wolter, B. Fuhrmeister, A. Reiners, G. Peres: Coronal explosion on CN Leo 




XMM-Newton OM-U 




5 10 15 20 25 30 

Seconds (since UT 23:46:25) 



Fig. 1. XMM-Newton EPIC raw X-ray (sum of PN,MOSl,MOS2 detectors) and OM U band light curves (bottom panels), UVES 
blue and red channel (upper panels) during the first 30 seconds of the impulsive flare on CN Leo ; all times are in UT seconds relative 
to May, 19 th , 2006, 23:46:25; the inserts show the same light curves covering 1000 seconds relative to May, 19 th , 2006, 23:39:40 in 
logarithmic scale including some quiescent emission and the full flare light curve; the light curve portions displayed enlarged are 
gray shaded. The X-ray light curves exclude events recorded in the center of the point response function to avoid pile-up effects; 
the total X-ray count rate is almost twice as large. The binning of all data is 1 sec, except that of the X-ray data in the insert, which 
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fraction of the optical light is directed into an exposuremeter for 
each of the two spectral bands ("blue" and "red"), thus provid- 
ing broad-band photometry with a time resolution of one second. 
These data are, however, taken mostly for engineering purposes, 
and are not corrected for background or airmass. For our anal- 
ysis we used the UVES pipeline products with the wavelength 
calibration carried out using Thorium- Argon spectra resulting in 
an accuracy of ~ 0.03 A in the blue arm and ~ 0.05 A in the 
red arm. Absolute flux calibration is based on the UVES master 
response curve and extinction files provided by ESO. 

3. Results 

3.1. X-ray and optical light curves 

In Fig. Q] we plot - from bottom to top - the recorded 
XMM-Newton X-ray light curve (i.e., the sum of the 
PN+MOS1+MOS2 detectors), the XMM-Newton Optical 
Monitor (OM) U band light curve, and the UVES blue and red 
light curves for part of our observing run. In each panel we show 
a linear light curve covering 30 seconds around the flare onset, 
while the respective inserts show the logarithmic light curves 
of the whole flare event including some quiescence. The most 
notable feature in these light curves is the enormous, simultane- 
ously occurring intensity increase in all considered energy bands 
starting at around UT = 23:46:20. The relative increase in terms 
of flare peak to pre-flare flux is > 100 in X-rays and « 750 in 
the UVES "blue" band. Even in flare stars flares with magnitude 
increases > 7 as shown in Fig. Q] are quite rare, and we are actu- 
ally not aware of any examples in the literature sampled with the 
time resolution and spectral coverage of our data. The temporal 
alignment of the light curves in different bands ought to be accu- 
rate to at least 0.5 seconds and hence below the bin width. Fig.Q] 
shows four remarkable properties of our multiwavelength light 
curves: First, the flare onset appears to be earlier in the optical 
wavebands, in particular, there is no evidence for any significant 
soft X-ray "precursor" activity, while there is some "pre-flash" 
visible in the U band; second, the rise from quiescence to flare 
peak in the optical is 750 fold and occurs within « 10 seconds; 
third, both the "red" and "blue" UVES light curves show two 
peaks separated by 5 seconds; and fourth, the first optical peak 
is accompanied by a soft X-ray peak. Note that this first soft X- 
ray peak is very sharp, lasting at most ~ 2 seconds, and does 
not coincide with the main soft X-ray flare peak, which in fact is 
quite flat and occurs about 200 seconds later (cf., insert in bot- 
tom panel of FigJT}. 



CN Leo flare onset 




-2TAPEC, red. x ! =0.89 
-powerlaw, red. x 2 -2.00 



Energy [ keV ] 



Fig. 2. EPIC pn spectrum of burst between UT 23:46:36 to 
UT 23:46:44 (black data points) together with model fit assum- 
ing a thermal spectrum (blue curve) and a power law (red curve). 
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Fig. 3. Comparison of burst spectrum between UT 23:46:36 to 
UT 23:46:44 (black data points) with postburst spectrum (red 
data points) between UT 23:46:44 to UT 23:47:44; a power law 
fit with photon index y — 2 is also shown. 



3.2. The spectral nature of the X-ray burst 

The crucial question concerning the X-ray light curve is the na- 
ture of the soft X-ray peak at time T = 15 sec in Fig[T] We em- 
phasize that this peak is statistically highly significant. We find a 
total of 105 counts in the two adjacent 1 sec bins around the opti- 
cal flare maximum. A smooth interpolation of the flare rise light 
curve leads us to a "background" estimate of w 55 counts in those 
two bins, implying an extremely high statistical significance of 
the 2 second flare peak. Accepting then the X-ray and optical 
peak and their close temporal association as real, on the one hand 
would suggest a non-thermal origin of the observed X-ray emis- 
sion, while on the other hand, non-thermal X-ray emission has 
never been conclusively demonstrated from stars (lOsten et all 
120071) . Also, impulsive soft X-ray ra diation is known fro m the 
Sun only on much longer time scales (Hudson et al., 1994). For 
the Sun, non-thermal, often impulsive X-ray emission is rou- 



tinely observed, but below ~ 25 keV the discrimination between 
thermal and no n-thermal emission becomes exceedingly difficult 
(iKahlerl fl975h . On the other hand, assuming a thermal origin of 
the first X-ray peak makes the observed fast decay time of ~ 2 
sec difficult to reconcile with the usually invoked conductive or 
radiative cooling time scales. 

Unfortunately the spectral analysis of the EPIC pulse height 
data in the respective time interval from UT 23:46:36 to 
UT 23:46:44 is plagued by the low SNR in this short stretch 
of data. The rebinned EPIC-pn spectrum of this time interval 
is shown in Fig. [2] together with a best fit assuming a thermal 
(blue curve) and a power law (red curve) spectral shape. For 
the thermal fit we assumed a two-temperature model with so- 
lar abundance, incidentally, a one-temperature component model 
with variable abundance also yields fits of comparable statisti- 
cal quality. The fit quality of a power law fit (x 2, d = 2.0) is 
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clearly far inferior to that of the thermal models (x 2 red = 0.9) , 
and a an improvement in quality of a power law fit can only be 
brought about by the introduction of substantial amounts of cold 
absorption (Nh ~ 4 x 10 21 cm -2 ), which we regard as utterly un- 
physical in the context of CN Leo. Therefore we clearly prefer a 
thermal interpretation of the "burst" spectrum purely on obser- 
vational grounds. This thermal interpretation is supported by the 
X-ray spectrum between UT 23:46:44 to UT 23:47:44, immedi- 
ately after the burst, displayed in Figf3] together with the power 
fit producing an acceptable fit to the burst spectrum in Figf2] 
Clearly, this spectrum cannot possibly be fit with a power law 
spectrum, the iron L complex and Ne lines in the energy range 
0.7-1.2 keV prevent a power law description; formally, a power 
law fit with variable slope results in (x 2 red — 3.32), while a ther- 
mal spectrum yields (y 2 prf = 1.01) and thus we conclude that the 
spectrum between between UT 23:46:44 to UT 23:47:44 is def- 
initely thermal. Going back to Fig.|2]we realize that the largest 
power law fit residuals are found near 1 keV, precisely in that 
spectral range where the thermal nature of the post-burst spec- 
trum becomes most apparent. If we assumed the burst spectrum 
to be non-thermal, the best fit power law slope of y — 2 of the 
burst spectrum (cf., Figf2]l would have to agree only coinciden- 
tally with the "slope" of the thermal spectrum in the post burst 
spectrum (cf., Fig|3). Therefore we conclude that the EPIC pulse 
height spectra provide no proof for a non-thermal character of 
the observed X-ray emission during any time in the flare. 




Fig. 4. VLT/UVES spectrum (between UT 23:30:35 to 
23:47:15) covering the wavelength range between 3300 - 3850 
A . Data points are shown as blue histogram, black body fits with 
temperatures of 1 1318 K in red and for 10 000 K and 15 000 K 
as dashed lines. Balmer line emission is visible at longer wave- 
lengths, but at shorter wavelengths most of the flux resides in 
continuum. 



3.3. Optical burst spectrum 

A high-resolution UVES spectrum covering the time interval 
UT 23:30:35 to UT 23:47:15 covering the wavelength range be- 
tween 3200 A - 3850 A has been taken (see FigEJl. While this 
spectrum covers also periods of quiescence, it is clearly domi- 
nated by the flare and can be taken as a representative mean spec- 
trum of the first minute of flare emission. Over this wavelength 
range the emission is almost linearly decreasing with wavelength 
without any major emission lines except the Balmer lines at the 
long wavelength end; the bulk of the chromospheric emission 
lines appears only in later spectra. The slope of the best fit black- 
body spectrum spectrum corresponds to an effective temperature 



of ~ 1 1 320 K; other fits with temperatures fixed at 10 000 K and 
15 000 K are also shown. Converting then the observed spectral 
normalization to an emission area (using the known distance to- 
wards CN Leo) results in an area estimate of su 4.4 x 10 18 cm 2 , 
while the other temperatures lead to values between 1.7-7.1 x 
10 18 cm 2 . Some caution is warranted in this estimate: As is ev- 
ident from Fig[l] significant photospheric temperature changes 
take place during the first 60 seconds, which cannot be resolved 
in the UVES spectrum. The derived temperature is therefore a 
mean and probably lower than the achieved maximum tempera- 
ture. Further, no calibration spectra of standard stars were taken 
during the same night and using the spectra of different standard 
stars taken during different nights leads to errors in the abso- 
lute calibration of up to 50 %. We therefore conclude that the 
emission area responsible for the observed blue optical emission 
was in the range 1-10 x 10 18 cm 2 . An independent estimate of 
the size of this emission area can be derived from the UVES 
photometer data. Using the mean temperature derived from the 
spectrum and the count rate during the spectral exposure, we 
can convert the observed count rate into an energy flux (from a 
comparison with other stars) and derive in this fashion an emis- 
sion area of w 8x 10 18 cm 2 , in good agreement with the estimate 
from the spectrum. We thus conclude that the photospheric flare 
peak temperature was in excess of 1 1 000 K and covered a flare 
area A// of 1 - 1 X 10 18 cm 2 and obviously one expects the X-ray 
emission to arise from a similarly-sized region. Using the pile-up 
corrected PN light curve we can use the observed flare count rate 
in the energy band 0.5 - 8 keV to compute an (isotropic) radiative 
loss of a 3.8 10 31 erg; we estimate that about two thirds of the 
overall radiative losses are contained in this energy band. During 
the first X-ray peak, the recorded X-ray luminosity is estimated 
to be ~ 5 x 10 28 erg s" 1 (Fig.©. 



4. Modeling and Interpretation 

4.1. Model ansatz 

To address the question whether our observations and the above 
derived physical parameters are consistent with the expected hy- 
drodynamic response of coronal plasma to a fast transient heat 
flux we performed a hydrodynamical simulation of an initially 
hydrostatic loop with prescribed length L — 4 x 10 9 cm and apex 
temperature T max = 2x 10 6 K subjected to impulsive heating. For 
our hydrod ynamical simulations we use the Pa lermo-Harvard 
(PH) code jPeres et al.L Il982t iBetta et al.L I1997I) . which solves 
the time-dependent 1-D hydrodynamic equations on an Eulerian 
grid. The energy equation includes the effects of optically thin 
radiative losses, heating and thermal conduction, treated in the 
Spitzer form, which should apply in a high-density medium, 
where the heat pulse is deposited, The heating term Q consists 
of two contributions, a spatially and temporally constant term 
to produce a quiescent corona, and a spatially and temporally 
varying term to describe a flare. Since the microphysics of the 
heating is not well understood, we use an ad hoc parameteriza- 
tion of the form Q/i - Hof(t)g(s). For g(s) we use a Gaussian 



g(s) = 



■.exp( ^1 ), centered at so 



2 x 10 s cm (as 



measured from the loop footpoint) and cr = 10 8 cm, for f(t) we 
use a symmetric triangular profile such that f(t) differs from 
for a time interval tj wa , and assumes a maximal value of unity 
and Hq = 7200 erg cirT 3 s _1 . The main difficulty for a numerical 
solution of the corona problem is the treatment of the transition 
region, where temperature and density gradients are extremely 
large. These numerical difficulties are aggravated by the fact 
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Length along lopp (1000 km) 



Fig. 5. Calculated run of temperature (upper panel), density 
(in units of 10 9 cirT 3 ; medium panel) and emission measure 
weighted cooling function (bottom panel) for the first 20 sec- 
onds of a coronal explosion. The initial state (a quiescent loop 
with T mflA = 3 x 10 6 K) is shown as dotted line, the evolution at 
times t = 2s, 4s, 8s and 10s, during which heating takes place, 
is shown as solid lines, that at times t = 13s, 16s, and 19s as 
long-dashed lines. 



that during the flare evolution the position of the transition re- 
gion changes its spatial location. The PH code uses a re-adaptive 
grid, which insures a maximal relative change of 10 % in all 
physical variables on adjacent grid points, provides a smoothly 
varying grid size over the whole computational domain and sam- 
ples the whole region of interest with typically 500 - 550 grid 
points. At the footpoints the model loop is anchored to a dense 
chromosphere, the s tructure of which follows the treatment by 
IVernazza et al.l dl98ll) and in which the energy balance is main- 
tained throughout the simulation. We assume the loop to be sym- 
metric with respect to its apex, so that the evolution is computed 
along one loop half and symmetry boundary conditions are im- 
posed at the apex. The lower boundary is maintained at a fixed 
temperature at the chromospheric temperature minimum. 

4.2. Explosion modeling 

We specifically consider a simulation with Td ep = 10 sec and a 
total deposited heat fluence E tot = 1.8 10 13 erg/cm 2 , assumed 
to be deposited at a height of 2000 km above the photosphere. 
The results of our simulation are shown in Fig. [5] where we plot 
the run of density and temperature vs. loop length for the initial 



state and every 2 seconds into the flare evolution. The main ef- 
fect of the transient heating is a rapid temperature increase of the 
"old" (i.e., pre-flare) corona and transition region of the loop and 
parts of the "old" chromosphere to temperatures of 20 - 40 MK. 
Because of the high-density of the latter, a high-pressure region, 
sandwiched between the low-pressure "old" corona and low- 
temperature unheated chromosphere/photosphere, is formed at 
the bottom of the "old" transition region. An expansion wave 
starts propagating upwards with speeds in excess of 1000 km/s 
to increase the density of the "old" corona ("loop filling"). This 
explosive filling of the loop occurs on a time scale of 10 seconds, 
comparable to the sound crossing time of the structure. Once the 
loop is rapidly filled, the material is more gently accumulated 
and compressionally heated, leading to a slower rise in emission 
measure and temperature. At the same time, the high-pressure 
"new" transition region moves inwards into the chromosphere 
with speeds of ~ 30 km/s; a narrow density spike separates in- 
ward moving outer regions from not yet accelerated inner re- 
gions. Once the heating stops, the temperature in the lower loop 
portions drops through conductive losses into the chromosphere, 
while the temperature at the loop top still increases through com- 
pressional heating. From the runs of temperature and density we 
compute the temperature-dependent emission measure and cool- 
ing in the energy range 0.5-8 keV, displayed in Fig. [5] Obviously 
during the first w 10 seconds the main contribution to the X-ray 
flux comes from the bottom parts of the loop, while at later times 
the evaporated plasma dominates the overall coronal emission. 



6x10 




5 10 15 20 
Seconds (since start of flare heating) 

Fig. 6. XMM-Newton X-ray light curve (in erg/s) overplotted 
over hydrodynamic model prediction; time is with respect to the 
modeled start of the energy input. Note the good agreement be- 
tween data and model although no best fit was performed. Also 
note that the event at t = 13 sec might represent a further soft 
X-ray burst. 



4.3. Comparison to observations and interpretation 

The theoretical curves displayed in FigfS] can be readily com- 
pared to our XMM-Newton data. It is straightforward to com- 
pute the total X-ray flux assuming some loop cross sectional area 
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and in Fig. [6] we juxtapose our model predictions (solid line) 
computed under the assumption of an emission area of A/; = 
10 19 cm 2 and the observed X-ray luminosities computed from 
the pile-up corrected PN lightcurve using a count-rate-to-flux 
conversion of 1.8 x 1CT 12 erg/cm 2 /count and the known distance 
towards CN Leo. Although no best fit (in terms of the temporal 
model evolution) has been attempted we note the good agree- 
ment between theoretical expectations and observations. The X- 
ray peak coincides with the end of the heating phase and the 
spike in the light curve at t = 13 sec may represent another soft 
X-ray burst. 
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4.5. Thick target bremsstrahlung 

As is clear from the above description the energy input in our 
modeling occurs in an ad hoc fashion by specifying an arbi- 
trary spatially and temporally varying heating function. Since 
the precise nature of the heating physics is unknown, this ap- 
proach seems justified. Of course, the physical picture we have 
in the back of our minds is that of energetic particles (elec- 
trons) penetrating deep into the chromosphere and upper photo- 
sphere, dissipating their energy there and heating up the adjacent 
atmosphere. Since the UVES light curves are well co-aligned 
with the initial X-ray light curve and since the former emission 
comes from relatively cool material far extending into the red 
wavelength band, we require fairly energetic electrons to ener- 
~ gize these photospheric layers. Unfortunately, the precise prop- 

- erties of these electrons and, in particular, their energy spectra 
J and spectral cutoffs are only poorly known for solar flares and 

essentially unknown in the case of stellar flares. At any rate, 
one expects thick-target bremsstrahlung from this non-thermal 

- particle population , but at what level ? To this end we follow 
: the ca lculations by Svratovskii a nd Shmeleval i 1972b and[Brown 
: d 19731) . who considered the bremsstrahlung X-ray spectra for 
~ bremsstrahlung emitted by a population of electrons injected into 
z the flare region with an energy distribution function of the form 
* 

1 F(E) = F ( A y . (1) 



20 40 60 80 

Seconds (since 85590) 

Fig. 7. XMM- Newton EPIC light curve (bottom panel) and 
mean photon energy (top panel) vs. time for the first 100 sec- 
onds around flare onset; all times are relative to May, 19 , 2006, 
23:46:20. The mean photon energy has been derived from all 
photons with registered energies > 500 eV; each data point rep- 
resents the mean energy of 200 photons. 



J The approach taken by iBrownl (119731) is analytical, the most 
recent calcula tions of thick-target bremsstrahlung are those by 
^Holmanl d2002l) . who computed microwave and X-ray spectra 
for such power law electron distributions using the appropri- 
ate relativistic cross sections and considerin g the effects o f low 
and high energy cutoffs. As demo nstrated bv lBrownl ( 1 19731) and 
Svratovsk ii and Shmeleval d!972l) . an electron power law distri- 
bution of the form given in eqQ]leads to a power law distribution 
of the observed photon flux density 1(e) of the form 



/(e) = /(»( 



(2) 



4.4. Comparison to observations and interpretation 

The above interpretation of the XMM-Newton light curve is sup- 
ported by the observed run of the evolution of count rate and 
mean photon energy for the first 100 seconds into the flare shown 
in Figj7] As is clear, the mean photon energy and therefore the 
mean X-ray temperature are decreasing, implying that the X- 
ray emitting plasma is cooling. In the downward direction con- 
ductive cooling is the dominant cooling process; computing the 
cooling time scale at the heat pulse center through conduction 
via T cond = jItt we find T cond * 1.7 sec using n = 7 x 10 12 
cm" 3 , L = 1 x 10 8 cm and T = 3.4 x 10 7 K. The rapid cooling 
in the lower section of the loop leads to a pressure deficit ac- 
companied by a density decrease and ensuing decrease in X-ray 
emission (cf., Fig|5]l, driven by the fast evaporation front mov- 
ing ahead. The short lasting soft X-ray burst shown in FigQ]is 
thus explained with a highly coherent plasma evolution, i.e., by 
flaring plasma confined in a single magnetic flux tube and by 
assuming that a major fraction or all of the observed energy re- 
lease of > 3.8 10 31 erg occurred on a rather short time scale of 
10 seconds or less. We therefore conclude that we are directly 
witnessing the process of evaporation as a result of the explosive 
heating of the lower chromosphere. 



i.e., again a power law distribution. The respective power law 
indices and normalizations are related through 



B = y-\ 



FqEq 



(3) 



AB(fi- 1,1/2)' 

B(J3- 1, 1/2) in eq.[3]denotes the Beta-function, and the constant 
A is given by 



A = - 



8 ri 



3 137 7re 4 Z 2 A' 



(4) 



where the quantities c,r t „m,e,Z and A denote speed of light, clas- 
sical electron radius, electron mass, electron charge, ion charge 
number, and Coulomb logarithm respectively. Finally, the total 
energy content E m of the electron population in the energy range 
from Eq to infinity can be computed from 



- - 



IqBE 



AB(B- 1,1/2) 



(5) 



Thus, as well known from solar physics, the power law index 
of the electron distribution as well as the normalization of the 
electron distribution and therefore its total energy content can be 
inferred from the observed photon power law index B and the 
observed photon fluence Iq, if the photons are indeed generated 
by the thick target bremsstrahlung mechanism. If we assume - 
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for argument's sake - a thick target bremsstrahlung origin for 
the observed burst, we use the best-fit photon spectral index of 
P ~ 2 and the burst integrated photon fluence at 1 keV of 4.8 
x 10 37 photons/keV to compute the total energy content E p i of 
the power law electron population of 2.4 xlO 35 ergs with a low 
energy cutoff of 1 keV; assuming a low energy cutoff of 0. 1 keV, 
we instead find E p i - 2.4 x 10 36 erg. Regardless of the assumed 
low energy cutoff, the thus derived total energies are a few orders 
of magnitude above the observed soft X-ray energy output of 3.8 
xlO 31 erg. In other words, under the assumption of non-thermal 
bremsstrahlung we would have to assume that unlike the solar 
case the bulk of the energy contained in the non-thermal elec- 
tron population goes into other (unobservable) forms of energy. 
Applying Occam's razor we thus conclude that both observa- 
tions and theory suggest that the observed soft X-ray emission 
is not due to thick-target bremsstrahlung. At the same time we 
have shown that thick-target bremsstrahlung, which is of course 
expected to be produced at some level, can be easily accommo- 
dated with in the existing observational framework, even if its 
low energy cutoff were to extend down to below 1 keV and there- 
fore our theoretical modeling is fully consistent with the existing 
data. 

4.6. Physical consistency of modeling 

Does our theoretical modeling capture all of the important 
physics of the flare explosion ? A basic assumption in our model- 
ing is the 1-D character of the problem. Enormous pressure dif- 
ferences exist between "bottom" and "top" of the flaring plasma, 
while the magnetic field imposes a strong 1 -dimensional direc- 
tionality of the plasma along the field lines. Inspecting the maxi- 
mally occurring gas pressures, we find values of up to 10000 dyn 
cm -2 , requiring a magnetic field of > 300 G for confinement. The 
mean photospheric field in CN Leo was si multaneously mea- 
sured to exceed 2000 G (Reiners et al., 2007), thus plasma con- 
finement is no problem and our 1-D modeling should be correct. 
The released energy is assumed to come from magnetic recon- 
nection. An efficient energy conversion requires the incoming 
flow velocity to be less than 10 % of the Alfven speed v^ and 
hence the time scale Tdt ss for an efficient magnetic field dissipa- 
tion on a spatial scale L becomes tj,,, ~ g-j^ . Using L ~ 2 x 10 9 

cm and t^ I!S ~ 5 seconds, we find v^ ~ 4 x 10 9 cm/s, which 
appears unrealistically large. However, since we expect coronal 
fields on the order of w 300 G and particle densities of as 5 10 8 
cm -3 , we find Alfven velocities of v^ ~ 3 x 10 9 cm/s, and there- 
fore the assumption of Alfven speeds in excess to 10 9 cm/s is 
reasonable. 



5. Discussion and conclusions 

Heating by so-call e d "microflares" and "nanoflares" 
dcEgill and Klimchul I200H iKlimchukL [2006) is a popu- 
lar hypothesis to explain the heating of the corona of the Sun 
and that of other stars. Their nanoflare model assumes that 
any coronal loop is composed of many unresolved strands of 
magnetic flux and that these strands are "heated impulsively 
by a small burst of energy (which, for convenience, are called 
a nanoflare, although the range of energies is completely 
arbitrary". Further, the strands are thermally isolated from each 
other and can be viewed as eleme ntal magnetic flux tubes with 
a small value for the plasma beta (Cargill and Klimchuk, 2004). 
An individual heating event on such a strand leads to a dynamic 
evolution of the affected plasma located on this particular 



strand. In a typical observation of a star and even in a typical 
observation of the solar corona, a multitude of heating events 
occur more or less simultaneously on different strands located in 
the field of view or are integrated over during a typical exposure 
time. Therefore in this picture any solar or stellar coronal 
observation represents a mean of the heating and cooling 
history of possibly a very large number of individual energy 
releases. For a computation of the global properties of solar and 
stellar coronae in terms of the distributions of plasma densities, 
temperatures and flow velocities, a knowledge and modeling of 
the relevant heating and cooling processes is mandatory. Such 
models specifically assume an impulsive heating of the strands 
of magnetic flux, which can be recognized only in data with 
very high time resolution (~ 1 sec), and then cool, first, b y 
conduction and then by radiation (Car gill and Klimchukl 12004). 
This is precisely what appears to have happened on CN Leo: 
An energy release on a time scale of less than 10 sec in a single 
loop has to our knowledge never been observed in any star nor 
on the Sun. Therefore, with our high time-resolution soft X-ray 
and optical observations of the impulsive flare on CN Leo - by 
its energetics definitely not a nanoflare - we may have isolated 
the relevant physics for nanoflare heating of the corona of the 
Sun and that of the stars. 
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